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A  computer  based/controlled  PIC  solid  state  electronic  system  is 
presented  in  its  entirety.  Its  operation  is  described  and  evaluated  both 
as  a  plasma  diagnostic  instrument  and  a  basic  reactor  power  measurement 
system.  The  system,  in  conjunction  with  two  sets  of  matched  ionization 
chambers  (one  with  a  93%  enriched    U  coating  the  other  without),  is 
evaluated  in  terms  of  response  in  radiation  fields  of  I0C  tc  3  x  10  R/hr 
at  temperatures  ranging  from  25°C  to  475°C.  Two  different  sets  of  cham- 
bers were  used,  one   containing  Ar-5%N0,  the  other  Ne.  Gamma  compensation 

to  within  %   -■•_  5%   proved  feasible,  at  c   qivin  temperature  from,  10'  to  3 

-,.,5  .--,. 
x  10  R/nr. 

"he  ensmbers1  responses  altered  over  the  temperature  r?j\nc   used. 

The  PIC  s^nai  for  neon  at  3  x  10  R/hr  zna   ?5°C  was  2.03V  where  at  432°C 

it  was  4,75V,  For  the  Ar  -  5%   M0  chambers  under  the  same  conditions 

the  signal  increased  from  1  .9V  to  ? .  74  vol  ts .  The  Ar-N,-,  chamber  temperature 

response  was  expected.  The  Ha   temperature  resoor.ee  was  constant  up  to 


iy. 


282°C  as  predicted.  However,  at  this  temperature  the  signal  began  to 
shew  a  slow  increasing  trend  which  continued  up  to  the  maximum  temperature 
used,  475°C.  This  requires  more  study. 


CHAPTER  I 
INTRODUCTION 


The  need  for  dependable  and  accurate  monitoring  of  neutron  flux 
(reactor  power)  presents  the  nuclear  engineer  with  cue  of  his  most  diffi- 


cult design  tasks..  For  safety  reasons,  commercial  reactors 


of  all  types 


require  the  flux  to  be  monitored  over  st  least  twelve  decades  tfith  +10% 
accuracy,  ''ins  range,  even  under   optimum  environmental  conditions, 
requires  a  minimum  of  two  detector  operational  modes.  However,  in  all 
cases,  three  nodes  are  used  to  insure  over! 20  (redundancy)  and  thus  the 
safe  operation  of  the  reactor.  In   addition  to  the  range  requirements, 
the  detectors  must  operate  consistently  and  accurately  'Ddr   the  jig s t 
adverse  environmental  conditions.  These,  depending  on  the  situation, 
include  high  temperatures,  excessive  baccgreund  radiation  levels,  and 
high  neutron  flux  levels,  ever'  extended  periods  of  time  (months). 

For  the  most  part,  adequate  monitoring  systems  exist  for   those 
nuclear  plants  which  are  in  us?  at  the  present  time   However,  the  i;mits 
of  these  current  systems  are  exceeded  when  one  attempts  tc  ^poly  them  to 
the  nee.d~-   of  the  fast  breeder  reactor  or  the  new  generation  of  large  core 
light  water  reactors.  These  now  plants  reu-^rk   fast  response,  in-core 
detection  systems.  Unfortunately,  the  in-core  environment  is  extreme  in 
both  the  temperature  and  flux  domains  and,  as  yet,  no  system  has  func- 
tioned adequately  under  such  conditions. 

To  delineate  the  exact  nature  of  the  problems  faced,  one  must  care- 
fully examine  the  reactor  power  measurement  systems  preser.tlv  beinq 


utilized.  Such  an  examination  will  not  only  suggest  where,  and  possibly 
how,  improvement  can  be  made  in  current  systems,  but  it  will  also  point 
to  the  pressing  need  for  a  full  range,  high  temperature  single  sensor 
system. 

The  first  division  between  neutron  sensing  systems  comes  from 
whether  they  utilize  in-core  or  ex-core  sensors.  Examining  figure  1-1  we 
see  why  such  a  division  exists.  There  is  a  difference  of  a  factor  of  10 
in  temperature  and  a  factor  of  10  in  gamma  and  neutron  fluxes  between 
the  two  regions.  It  is  apparent  that  those  systems  designed  for  in-core 
monitoring  must  be  an  extremely  hearty  breed.  For  this  reason,  only  one 
of  the  four  major  reactor  types  (PWR,  HTGR,  LMFBR,  and  BWR)  has  the 
neutron  flux  level  safety  control  monitors  in-core.  The  BWR,  due  to  its 
large  volume  xnd   particular  design  characteristics  (large  cruciform 
control  rods  and  boiling  water)  requires  in  core  flux  monitoring  for 
safety  reasons.  This  requirement,  coupled  with  the  "lack  of  dependable 
in -core  systems,  together,  comprise  one  of  the  most  difficult  problems 
the  BWR  engineer  faces  today. 

Figure  1-1  also  serves  to  define  the  upper  limit  operating  specifi- 
cations for  both  the  in-core  and  ex-core  environment.  The  lower  limits 
corresponding  to  reactor  shutdown  are  depicted  in  figure  1-1  as  well. 
The  great,  range  in  reactor  power  makes  the  use  of  a  single  sensor  and 
circuit  impossible  with  current  technology.  Thus  a  detector  of  any  given 
design  is  of  use  only  over  a  specific  part  of  the  flux  range  -md   must  be 
complimented  with  other  sensor-circuit  designs.  Overlapping  one  or  two 
decades  of  each  design  sacrifices  part  of  their  useful  range  but  insures 
smooth  transfer  of  control  and  safety  functions  from  one  control  circuit 
to  the  next. 
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Figure  '1-1.  Typical  environmental  profile  for  neutron  sensors  in  a  PWR 
at  full  power  ana  at.  shutdown 


Three  ranges  are  used  to  divide  neutron  flux  levels  found  in  reactors 
from  startup  to  full  power:-  the  source,  intermediate,  a^.d   full  or  operat- 
ing power  range.  Figures  1-2  and   1-3  show  the  limits  of  these  ranges, 
along  with  the  typical  neutron  detectors  used  in  common  ex-core  and  BWR 
in- core  control  systems. 

It  is  apparent  from  figure  1-3  that  the  source  range  monitoring 
system  must  Jtilize  sensitive  detection  methods.  Proportional  and  fission 
counting  systems  offer  the  greatest  sensitivity  with  maximum  gamma  dis- 
crimination. Three  neutron  sensitive  materials  are  commonly  used  in  the 

10        3         10 
proportional  counters;   BF-,  gas,  He  gas  and   B  as  a  'lining.  Each  has 

res  advantages  and   disadvantages.  '  BF,  offers  high  sensitivity  but  the 

necessary  high  polarizing  voltages  cause  rapid  degradation  of  the  gas  in 

3 
high  flux  environments.  The  He  gas  proportional  counter  offers  greater 

sensitivity  and  stability  than  the   BF~.  but  has  a  smaller  Q  value, 

10 
which  makes  the  gamma  sensitivity  proportionately  greater.  The   B  lined 

chamber  is  less  sensitive  than  the   BF-  but  more  stable  in  high  inten- 
sity fields.  Fission  chambers  have  their  electrodes  coated  with  uranium 

.  •   ;  ,  ■   .    A        •  235,  , 

highly  enriched  in   U. 

In  both  in~core  and  ex-core  systems,  the  minimum  allowable  count 
rate  for  safety  reasons  is  1  to  10  counts  ver   second.  Thus,  the  detec- 
tor's sensitivity  in  each  case  must  be  adjusted  to  insure  that  the  shut- 
down reactor  neutron  flux  results  in  a  count  rate  of  greater  than  1  count 
per  second.  Due  to  the  resolution  limitations  of  such  counting  systems, 
their  fastest  possible  response  is  on  the  order  of  10  counts  per  second, 
lo  the  point,  a  state  of  the  art  fission  counter  can  accurately  indicate 
neutron  flux  levels  over  six  decades  while  immersed  in  a  gamma  flux  as 
high  as  10  R/hr.  Although  counting  systems  hiv;j  their  problems,  no 
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Figure   >-2.     Typical   ranges  ami  detectors  used  in  ex-core  systems  to 
cover  the  source,   intermediate  and  power  ranges. 
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"igure  1-3.  Typical  ranges  and  detectors  used  in  in-core  systems  to 
cover'  the  source,  intermediate  and  power  ranges. 


other  system  yet  developed  offers  the  required  sensitivity  in  conjunction 
with  the  necessary  gamma  discrimination  that  they  do.  Counting  systems, 
when  stretched  to  their  limits,  cove?'  only  half  of   the  twelve  decade 
range. 

fhe  intermediate  range  instrumentation  covers  most  of  the  remaining 
six  decodes,  overlapping  two  to  four  decades  of  the  source  range  and  pert 
or  all  of  the  power  range.  The  sensors  used  to  spun  this  range  are  ion- 
ization chambers  of  various  designs  with  boron  or  fissile  electrode  coat- 
ings. When  simple  fission  (or  boron)  chambers  are  used,  the  associated 
signal  processing  circuitry  is  designed  so  that  the  resulting  system's 
output  is  some  measure  of  the  mean  square  voltage  (mean  of  the  squares  of 
the  deviation  from  the  mean),  abbreviated  MSV.  Since  Poisson  statistics 
govern  the  pulses  from  a  nuclear  radiation  sensor,  a  measure  of  the  mean 
square  voltage  is  a  direct  measure  of  the  mean  counting  rate.  This 
method  has  just  recently  been  developed  and  offers  at  least  three  impor- 
tant advantages  over  the  conventional  compensated  ionization  chambers:" 
increased  gamma  discrimination  (100  times  more  than  the  CIC) ,  improved 
operation  when  chambers  and  cables  are  exposed  to  elevated  temperatures, 
?in<i   more  efficient  use  of  sensors. 

Compensated  ion  chambers  (CIC)  are  the  sensors  commonly  used  to 
cover  the  intermediate  range.  Such  chambers  provide  signals  which  have 
the  gamma  induced  component  reduced  by  a  factor  of  20  to  100  times  ever 
that  of  a  conventional  ionization  chamber.  The  methodology  employed  to 
achieve  this  is  simple.  The  CIC  is  constructed  with  two  separate  sensi- 
tive volumes;  one,  having  the  confining  electrode  surfaces  coated  with 
neutron  sensitive  material,  is  sensitive  to  both  gammas  and  neutrons,  and 
the  other,  having  no  such  coatings,  is  sensitive  to  only  gammas.  The 


sensitivity  of  the  two  volumes  is  adjusted  so  that  they  ire   equally 
responsive  to  gammas.  Thus,  when  the  radiation  induced  ionization  from 
both  volumes  is  collected  and  subtracted,  the  resultant  signal  is,  in 
theory,  that  induced  solely  by  the  neutron  interactions.  Such  a  procedure 
is  necessary  because,  while  the  prompt  gamma  flux  is  proportional  to  the 
neutron  flux,  the  gamma  flux  due  to  radioactive  decay  is  not.  Thus,  all 
gamma  response  must  be  negated  as  much  as  possible.  This,  however,  is 
true  only  for  the  intermediate  range.  At  high  power  levels,  when  the 
neutron  field  is  much  more  intense  than  the  background  gamma  field,  no 
compensation  is  necessary.  Because  of  this  fact,  conventional  ionization 
chambers  are  used  as  the  control  monitors  from  1  to  100%  of  reactor 
power . 

The  use  of  a  gamma-compensated  detector  extends  the  reactor  control 
range,  compared  to  that  of  an  uncompensated  chamber,  by  approximately  two 
decades.  The  reason  this  extension  is  ;>o  small  lies  in  the  fact  that 
with  fixed  voltages  compensation  is  exact  at  only  one  given  reactor 
power-  level.  This  point  has  also  given  rise  to  the  recent  practice  of 
operating  with  fixed  voltages  end  designing  safety  systems  that  avoid 
total  control  dependence  on  more  than  two  deccdes  of  compensation. 

The  upper  limit  of  any  ionization  chamber's  operation  is  fixed  by 
either  recombination  affects,  which  cause  nonlinear  responses,  or  by  the 
inability  to  apply  sufficient  collection  voltage  to  the  electrodes.  The 
flux  level  at  which  either  of  these  occur  uepsnds  (^   the  overall  design 
of  the  sensur.  Leakage  current  through  the  insulatorss  due  to  the  applied 
voltage,  becomes  the  limiter  for  measuring  low  neutron  flux  levels.  In 
addition  to  insulation  leakaae,  the  lower  limit  of  the  ionization  chamber 


operation  may  be  affected  by  background  current  caused  by  the  activation 
of  chamber  materials,  neutron  sensitive  material  reaction  product  activity, 
or  bac k g re un d  alpha  current  for  f i s s  i o n  c h a mb e r s . 

It  was  mentioned  at  the  beginning  that  present  monitoring  systems 
are  adequate,  and  they  are,  but  just  barely  so.  Points  y;hicb  are  diffi- 
cult to  get  across  in  a  summary  of  this  kind  are  the  complexity,  failure 
rate,  and  field  engineering  problems  that  these  systems  present.  The  in- 
core  BKR  system  is  a  case  in  point.  How  is  one  to  calibrate,  normalize, 
and  keep  in  operating  condition  the  hundreds  of  detectors  present  in  the 
core  of  an  operating  reactor?  The  need  for  a  distinct  system  to  cover 
each  of  the  three  flux  ranges  makes  things  just  that  much  more  difficult. 

It  is  apparent  that  a  single  sensor,  full  range  reactor  power  measure- 
ment system,  which  could  cover  the  entire  flux  range  without  being  stretched 
to  its  limit,  would  vastly  improve  the  situation.  Such  a  system  should 
function  in  the  source  range  as  a  counter,  with  inherent  gamma  compensa- 
tion through  pulse  height  discrimination,  and  should  smoothly  switch  to 
the  intermediate  range  measurements  with  substantial  overlap  to  insure  a 
linear  transition..  At  intermediate  flux  levels,  the  system  should  be 
capable  of  both  gamma  conpens? cion  and  minimization  3f  leakage  current 
contributions  to  the  signal,  especially  at  the  elevated  temperatures 
encountered  in  the  in-core  environment.  Build  up  of  neutron  sensitive 
material  reaction  products  should  not  adversely  affect  the  system.  For 
power  measurements,  the  range  of  operation  should  not  be.   restricted  by 
recombination  effects.  To  be  of  use,  such  a  system  must  provide  a  linear 
(or  log)  output  over  the  entire  power  range  under   all  expected  environ- 
mental conditions.  The  Pulsed  Ionization  Chamber  (PJC)  technique,  which 
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is  to  be  the  topic  of  this  dissertation  research  report,  and  which  was 
developed  and  investigated  by  Ellis  and  his  students  at  the  University  of 
Florida,  appears  to  have  most,  if  not  all,  of  the  above  mentioned 
characteristics. 

The  promise  of  the  PIC  methodology  for  meeting  the  real  needs  and 
future  requirements,  alluded  to  above,  is  considered  to  constitute  more 
than  adequate  justification  for  undertaking  a  research  program  for  further 
developing  the  PIC  system  towards  practical  nuclear  power  reactor  applica- 
tions. Therefore,  to  better  impliment  the  concepts  and  demonstrate  the 
desired  operational  characteristics,  the  development  of  a  much  more- 
sophisticated  and  practical  solid  state  PIC  pulsed  high  voltage  and 
control  system  would  need   to  be  undertaken. 

Thus,  the  main  goal  in  initiating  the  research  described  in  the  sub- 
sequent chapters  of  this  dissertation  was  to  develop  and  evaluate  a 
single  sensor  compensated  PIC  system,  capable  of  full  range  in-core 
reactor  power'  measurements.  In  order  to  establish  an  initial  base  en 
which  to  develop  this  research  topic,  a  review  of  the  PIC  system's  basic 
operational  characteristics  and  previous  research  results  therewith  is 
presented  in  the  following  chapter. 


CHAPTER  II 
PULSED  IONIZATION  CHAMBER  OPERATIONAL  THEORY 


rhe  Pulsed  Ionization  Chamber  (PIC)  mode,  which  was  originally 
developed  for  plasma  diagnostic  purposes,  is  a  new  mode  of  operating  gas 
filled  ionization  chambers."  In  "its  initial  stages,  the  PIC  mode  was 
applied  to  the  measuring  of  ionization  densities  and  recombination  para- 
meters in  gas  filled  chambers  exposed  to  neutron  radiation.  The  logical 
deduction  made  at  that  time  was  that  if  the  PIC  output  was  known  as  a 

function  of  the  neutron  flux,  then  such  output  could  be  used  to  measure 

4 
unkown  fluxes.  This  was  proven  to  some  extent  as  reported  by  Markwell 

who  demonstrated  the  basic  PIC  performance  over  eight  decades  of  reactor 
power,  figure  2-1 . 

The  PIC  methodology  involves  the  periodic  application  of  a  single 
polarity  voltage  collection  potential  across  the  electrodes  of  an  ioniza- 
tion chamber.  Sufficient  time  is  allowed  between  high  voltage  pulses  for 
the  ionization  density  in  the  gas  filled  gap,  between  the  chamber  elec- 
trodes, to  approach  its  asymptotic  steady  state  limit.  The  application 
of  the  collection  voltage  results  in  the  collection  of  first  the  electrons 
and  subsequently  the  ions  of  the  equilibrium  ionization  density,  n.  The 
fact  that  the  electrons  arc  collected  approximately  a  thousand  times 
faster1  than  the  ions  makes  it  feasible  to  use  the  collection  of  only  one 
of  these  charged  particle  types  in  measuring  the  steady  state  ionization 
density.  Which  of  these  is  used  depends  almost  entirely  on  the  chamber 
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Figure  2-1.  Results  of  PIC  system  with  1  atm.  He  chamber. 
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design,  and  the  fill  gas  composition  and  pressure.  For  the  research 
reported  here,  the  ions  were  the  species  used  in  order  that  some  of  the 
problems  associated  with  electron  kinetics  could  be  avoided.  The  theory 
derived  and  presented  in  the  following  paragraphs  will  revolve  around  ion 
collection  rather  than  electron  collection  used  in  earlier  works. 

There  are  two  factors  basic  to  the  feasibility  end  usefulness  of  the 
PIC  system  for  the  field  of  radiation  measurement.  The  first  of  these  is 
that  the  ionization  density  in  a  ehan;b°r:,  when  it  is  exposed  to  ionizing 
radiation,  prows  rapidly  (10  to  100  milliseconds)  to  an  asymptotic  limit 
dependent  on  the  source  intensity.  The  second  factor  is  that  it  must  be 
possible  to  theoretically  relate  the  measured  voltage  signal,  due  to  the 
collection  of  ions,  to  the  asymptotic  ionization  density,  and  thus  the 
source  intensity. 

The  description  of  the  positive  ion  arid  electron  densities  growth  to 
their  asyrntotic  value  lies  in  the  following  gas  kinetics  equations ; 


an 
W 


<o..««. 


(2-1) 


3n._ 

-~-  =  v.  (D..vn. )  -  ex..  n  n.  +  S, 

where:  n  -  positive  ion  dens  icy  (cm  )s 
n  =  negative  ion  dens  icy  (cm  ), 
a,  =  first  order  electron  loss  coefficient, 
au  =  second  order  loss  coefficient, 
S  =  ionization  source  rate  (ion- pairs/second) , 
and    D  =  effective  diffusion  coefficient  which  is  functionally 
dependant  or,  the  radial  position  of  n  in  cylindrical 
geometry  because  of  the  variance  of  n  in  that  direction. 


(2-2) 


It  should  be  noted  that  these  coupled  equations,  when  used,  generally 
constitute  a  nonlinear  pair  of  equations  with  no  possible  analytical 
solution.  However,  for  the  densities  of  interest  in  this  study,  the 
equations  may  be  markedly  simplified.  The  two  ion  density  regions  o^ 
interest  are  those  controlled  by  the  free  diffusion  regime  and  the  ambi- 
polar-volume  recombination  regime. 

Second  order  loss  mechanisms  can  be  neglected  at  ion  densities  where 
fre.s  diffusion  dominates  and  the  asymtotic  form  of  equation  2-2  reduces 
to; 


M  n.  +  S  =  0 


(2-3) 


When  this  equation  is  applied  to  cylindrical  geometry  and  is  subject  to 
the  boundary  conditions, 


+      +■ 


where  a  and  b  are  the  outside  boundary  of  the  inside  electrode  and  the 
inside  boundary  of  the  outer  electrode,  respectively,  the  exact  solution 


n+(r)  -  4D~TnTb7a7  E(^-a2)^  "  ln  (-'^^r'   +  ^lr:b  ~  b  lna].  (2-4) 


At  ion  densities  greater  than  ^10  cm"  ambipolar  diffusion  is  con- 
trolling the  mechanism  for  spatial  distribution  of  charged  particles  and 
volume  recombination  becomes  the  major  loss  mechanism.  Because  of  ambi- 
polar diffusion,  the  ion-electron  diffusional  losses  become  approximately 


equal,  thus  validating  the  approximation  n,  =  n_  at  such  densities.  This 
approximation  leads  to, 

D  Vn.  -  »„nj  +  S  ■  0,  (2-5) 

which  is  the  decoupled  ion  density  equation  for  n>jQ  .  This  equation, 

when  written  in  cylindrical  coordinates,  is  the  Emden  differential  equa- 

g 
tion  which  cars  be  solved  using  only  numerical  methods . 

The  PIC  voltage  signal  amplitude  v(t  ),  at  the  time  t„,  when  all  of 

3    -1  r  c  c 

the  ions  are  collected,  for  a  large  RC  cathode  circuit  time  constant 
(large  with  respect  to  t  ),  and  for  coaxial  detector  electrode  geometry, 
is  analytically  and  experimentally  related  to  the  steady  state  ionization 
density,  n+,  in  the  chamber  gas.  by, 

rb 

v(tj=^/    rn.(r)dr,  (2-6) 

where;    1  =  length  of  the  chamber, 

C  =  cathode  circuit  capacitance, 
e  =  the  unit  of  electron  charge. 

The  final  link  in  relating  the  measured  vol  Lace  peak  signal  v(t  ) 
to  the  neutron  and  gamma  flux,  is  provided  by  the  equations; 


EaN  & 

-_*-  ,  (2-7) 


ana 


5  «  qly,  (2-8) 
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where;   E  -  the  average  energy  deposited  in  the  gas  per   neutron 
interaction, 

w  =  the  mean  eneryy  required  to  create  an  electron-ion  pair, 

a   ~  average  neutron  interaction  cress  section, 

iy  =  atom  density  of  the  neutron  sensitive  material, 
a 

$  -   the  flux  density  in  the  vicinity  of  the  detector, 

g  =  gas  ionization  efficiency  which  is  for  a  given 
radiation  field  arid  chamber  design 

I  :-  gamma  source  intensity. 

The  ion  source  rate,  S,  is  directly  relatable  to  the  relatively  easily 

measured  experimental  valve  I  ,  the  steady  state  ionization  current. 

1  c  is  measured  by  applying  a  constant  collection  potential  to  chamber 

electrodes  erA   measuring  the  current  between  those  electrodes  which 

result  from  the  ionizing  radiation. 

In  Darticular, 


eli  " 


(2-9) 


where,    e  =  the  unit  of  electronic  charges 

U  =  the  sensitive  volume  of  the  chamber. 
Equations  2-4  through  2-9  thus  provide  a  direct,  relatively  simple  rela- 
tionship between  the  PIC  output  signal  v(t  )  and  the  ionizing  radiation 
field  strengths,  and,  as  such,  serve  as  the  theoretical  basis  for  this 
entire  endeavor.  With  the  aid  of  the  author,  detailed  numerical  caicula 
tior.s  were  performed  by  Meravi  for  a  miniature  fission  chamber  filled 
with  neon  and  argon,  Some  of  the  results,  are  presented  in  figures  2-2 
through  2-5..  The  curves  in  figures  2-2  and  2-3  represent  the  radial 
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charge  distribution  in  an  argon   filled  chamber  for  low. and  high  ioniza- 
tion source  rates,  respectively.  The  severe  flattening  of  the  distribu- 
tion of  n  at  the  higher  source  rate  is  clue  to  volume  recombination, 
figures  2-4  and  2-5  depict  the  theoretically  determined  v(t  }  versus  S 
for  argon  and  neon  gas  filled  miniature  fission,  chambers.  In  both  cases 
the  response  is  at  first,  linearly  controlled  by  free  diffusional  losses. 
At  approximately  10"  ion  pair/sec,  ambipolar  diffusion  becomes  the  con- 
trolling factor  and  there  is  a  lateral  shift  in  the  linear  response  due 
to  the  more  rapid  ambipolar  diffusional  loss  mechanism.  Above  108  ion 
pairs/sec  volume  recombination  becomes  the  controlling  loss  mechanism  and 
the  response  becomes  second  order.  These  results  were,  to  some  extent, 
experimentally  validated  by  Markwell .   Figure  2-1  shows  measurements 
taken  in  the  University  of  Florida  Training  Reactor  (UFTR)  with  a  PIC 
system.  Over  seven  decades  of  reactor  flux  were  measured.  Both  first 
and  second  order  gas  kinetic  response  region  are  observed,  as  was  predicted, 
A  milliwatt  was  the  minimum  measurable  reactor  power,  due  to  the  noise 
and  time  jitter  of  that  PIC  system   The  upper  limit  of  10  Kw  was  then 
the  maximum  output  of  the  UFTR.  Note  that  care  must  be  taken  in  directly 
relating  the  experimental  and  numerical  data  shown,  since  the  former 
deals  with  the  collection  of  electrons  and  the  later  with  the  collection 
of  ions. 

The  effective  replication  of  the  characteristics  of  the  sets  of  data, 
resulting  from  either  gamma  or  neutron  caused  ionization,  indicates  the 
basic  dependence  of  the  PIC  operation  on  the  ionized  gas  kinetics,  regard- 
less of  which  radiation  induced  the  ionization.  Thus,  although  a  neutron 
sensitive  chamber  results  in  a  signal  which  is  comprised  of  both  neutron 


and  gamma  induced  ionization,  by  matching  this  Detector  with  one  that  is 
sensitive  to  gammas  only,  compensation  is  possible,  at  least  at  a  given 
temperature,,  as  was  shown  by  Cooper.   Some  of  his  experimental  results 
are  shown  in  figure  2-6. 

Conventional  in-core  ionization  chambers  suffer  greatly  from  high 
temperature  induced  leakage  current  effects  as  shewn  in  figure  2-7. 
The  departure  from  linearity  of  the  ionization  current-  as  a  function  of 
reactor  power  at  low  flux  levels,  results  from  temperature  enhanced 

leakage  current.  The  PIC  method's  ability  to  greatly  reduce  such  current 

4 
effects  was  shown  by  Ellis  and  Irnani  and  reported  in  Markwell's  thesis. 

Figure  2-8  shows  their  results.  The  change  in  signal  output  for  increased 

temperatures  was  shown  to  be  due  to  the  change  of  the  recombination 

coefficient  with  temperature.  However,  such  a  change  would  complicate 

the  application  of  the  PIC  methodology  to  reactor  power  measurements. 

p 
Fortunately,  however,  Sanders  results  indicate  that  this  drawback  could 

possibly  be  averted  in  chambers  using  neon  as  a  fill  gas,  since  neon's 

recombination  coefficient  was  reported  to  be  independent  of  temperature 

over  the  range  of  70°F  (25°C)  to  572^F  (300°C). 

Even  though  the  previously  developad  PIC  systems  were  adequate  for 

the  proof-of-principle  application  described  above,  they  fell  far  short 

of  being  prototypic  of  a  practical  reactor  power  measurement  system.  The 

9 
bases  for  such  a  practical  design  was  set  forth  by  Lllis.   The  principle 

technique  of  his  system  involved  direct  analog  signs'!  processing  in  com- 
bination with  direct  logic  system  gating  resulting  in  live  time  feedback 
range  and  operational  mode  control.  The  major  advantage  of  this  approach, 
in  terms  of  reactor  safety  and  control,  would  be  its  fast  response  time 
characteristics.  However,  it  would  suffer1  from  the  inherent  disadvantage 
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:igure  2-6.  Gamma  compensation  as  a  function  of  gamma  field  strength, 
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of  voltage  level  instabilities  and  drift  characteristics  of  analog  systems 
and  the  inflexibility  which  is  associated  with  hard-wired  signal  process- 
ing and  control  systems. 

State-of-the-art  digital  data  processing  avci   control  systems  are 
capable  of  speeds  more  than  adequate  for  reactor  application  with  the  PIC 
system.  Use  of  such  a  system  would  result  in  both  greater  system  stability 
and  adaptive  flexibility.  However,  such  a  high  speed  computer  for  develop- 
ing a  prototypic  system  might  incur  an  unnecessary  expense,  since  the 
performance  capabilities  of  a  high  speed  digital  computer  PIC  reactor 
control  system  could  be  adequately  demonstrated  through  the  use  of  a  more 
readily  available  digital  computer  having  less  stringent  response  and 
processing  characteristics,  Thus,  for  this  reason,  a  more  moderately 
priced,  but  adequate  unit  was  sought.  The  HP9821 ,  having  the  desired 
characteristics,  was  adapted  as  a  component  of  the  prototypic  digital 
computer  based  PIC  power  reactor  control  system  (PDCPIC)  which  was  de- 
signed and  built  for  this  dissertation. 

In  the  following  chapter  the  general  design  and  operational  charac- 
teristics of  the  PDCPIC  are  presented  in  detail.  This  is  roll  owed  in 
Chapter  IV  by  a  presentation  of  the  experimental  techniques  used  to 
evaluate  the  system  in  terms  of  gamma  compensation  and  detector  tempera- 
ture response:  two  important  operational  characteristics  which  needed 
further  study.  The  result?,  of  these  tests  are.-  then  evaluated.  Finally, 
in  Chapter  V,  the  author's  conclusions  are  put  forward  and  the  areas 
requiring  additional  research  are  stated. 


CHAPTER  III 
PIC  POWER  REACTOR  MEASUREMENT  SYSTEM" 


The  potential  benefit  of  the  research,  If  not  already  clear,  will  be 
clarified  in  the  following  general  description  of  a  PIC  full  range,  in- 
core,  single  sensor,  reactor  power  measurement  system, 

Three  operational  modes  are  to  be  used  to  cover  the  expected  power 
range:  count  rate,  compensated  PIC,  and  gamma  only  PIC  operation.  Based 
on  available  performance  data,  a  decade  or  more  overlap  between  these 
modes  for  safety  purposes  should  be  easily  attainable. 

The  count  rate  mode,  as  summarized  earlier,  is  the  best  low  flux 
measuring,  high  gamma  discrimination  method  available.  Integration  of  it 
into  the  PIC  scheme  requires  only  the  coupling  of  the  count  rafe  system 
to  the  detector  cathode  and  applying  a  constant  collection  bias  to  that 
electrode  as  depicted  in  figure  3-1.  Such  a  cor. figuration  would  cover 
the  in-core  flux  range  of  10°  to  10  neutrons/cm  /sec. 

The  CIC  mode,  as  previously  mentioned,  nas  some  serious  drawbacks; 
the  PIC  compensated  mode,  because  of  its  irsensitivity  to  leakage  current, 

should  offer  a  wider',  more  stable  range  of  response.  The  in-core  compen- 

7     13  2 

sated  PIC  system  would  cove"  a  6  decade  range;  !0  to  iO   neutrons/cm  / 

sec.  The  operation  of  this  mode  could  prove  to  be  quice  complex  and  will 

be  discussed  later. 

9  2 

Above  approximately  !0  neutrons/cm  /sec,  the  field  strength  from 

the  prompt  fission  gammas  begins  to  exceed  that,  resulting  from  radioactive 
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1  ]         2 
decay  of  the  fission  product  inventory.  Thus,  above  10  '  neutrons/cm  / 

sec,  the  output  of  the  gamma-only  section  of  the  compensated  chamber  can 
be  used  to  determine  the  power  level  of  3  reactor.  Depending  on  the 
overall  detector  design,  fluxes  as  high  as  10i;  could  be  measured  using 
this  technique.  Thus,  a  single  sensor  system  utilizing  the  modes  des- 
cribed, with  the  proper  sensor.,  would  indicate  reactor  power  over  the 
full  expected  range  of  operation,  encompassing  a  possible  neutron  flux 
range  of  10  to  10   neutrons/cm  /sec. 

To  design  and  assemble  such  a  system  from  the  basis  which  existed 
at.  the  initiation  of  this  program  would  not  only  have  been  difficult  but 
also  impractical.  While  verification  of  some  of  the  basic  PIC  responses 
were  performed,  the  entire  range  and  depth  of  its  response  had  not  been 
proven  as  a  whole.  For  instance,  neither  full  range  gamma  compensation 
nor  the  chamber's  response  to  temperature  variations  of  greater  than 
approximately  280° C  had  beer,  examined  fully. 

For  this  reason  the  prototypic  system  described  below  was  designed 
not  only  for  demonstrating  the  operational  characteristics  of  the  PIC 
instrumentation,  but  also  to  el  low  the  gamma  compensation  capability  to 
be  better  evaluated  and  the  temperature  characteristics  measurements 
to  be  extended  and  reverified.  The  assembling  of  such  a  prototypic 
system  not  only  indicates  the  feasibility  of  constructing  the  entire 
reactor  grade  system,  but  also  offers  a  vehicle  for  examining,  in  depth, 
all  the  PIC  methodology  characteristics, 

PIC  Electronic  System 
The  PIC  electronic  system  is  composed  of  5  basic  units;  the  high 
voltage  pulser  (HVP),  the  raw  data  acquisition  system  (PDA),  a  computer 


controlled  PIC  to  steady  state  mode  switching  system,  the  computer  inter- 
face modules  and  finally  the  HP9821  computer  system.  A  system  block 
diagram  is  given  in  figure  3-2,  Photographs  of  it  assembled  are  shown  in 
Figure  3-3.  With  the  exception  of  two  analog  to  digital  converters,  a 
picoaraneter,  certain  power  supplies,,  and  the  HP9821  system,  the  electronics 
were  designed  and  built  by  the  experimenter.  Because  of  the  newness  of 
the  design  approach  taken  and  the  resulting  improvement  over  all  previous 
PIC  electronics,  its  design  and  operation  are  carefully  outlined  below. 

The  solid  sUte  high  voltage  purser  is  depicted  in  figure  3-4.  There 
were  two  basic  design  goals  for  this-  unit.  The  first  was  to  apply  a  350 
volt  potential  -"cross  a  detector -cable  system  totaling  approximately  5C0 
pf  capacitance  in  approximately  500  nsec.  The  second  goal  was  to  return 
the  cathode,  which  was  essentially  grounded  in  order  to  meet  the  first 
requirement,  to  a  high  impedance  state  within  100  nsec  after  full  voltage 
applications  so  that  little  of  the  collected  charge  vnulc   be  lost.  Tnese 

two  goals  were  to  be  accomplished  using  solid  state  electronics.  Time 

.....  -     ,3,4,6.7,8,9    .  Ai  j 

.jittering  components  previously  used  .        such  as  mercury  wetted 

relays  «ere  to  be  avoided  at  all  cost. 

figure  3-5  gives  a  simplified  description  of  the  h:,sic  pulser  opera- 
tion. A  close  study  of  it  makes  following  the  circuit  description  below 
easier. 

The  integrated  circuit  (IC)  use-]  throughout  the  entire  system  is 
an  SN74121,  monostable  multivibrator.  This  14  pin  IC  is  triggered  by  a 
+5   to  0V  transition  input  at  pins  3  and/or  4.  Pins  I  and  6  provide  the 
resulting  negative  going  (+5  to  0V )  wd   positive  goir.y  i.u  to  5V)  output 
pulses,  respectively.  The  pulse  -iuths  are   controlled  by  the  capacitive 
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Figure  3-2.  Block  diagram  of  the  prototypic  digital  computer  based  PH 
reactor  power  control  system. 
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B.   SOME  OF  THE  INTERNAL  CIRCUITRY  OF  THE  MVP  AND  RDA, 
Figure  3-3.  Photographs  depicting  the  PIC  system. 
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>IC  equivalent  circuit'  with  pulsing  sequence 


lose,  across  pins  10  and  11  and  resistive  load  between  pin  9  and  >5  volt 
supply.,  Note  that  since  both  outputs  1  and  6  appear  simultaneously, 
utilizing  the  positive  rather  than  the  negative  output  to  trigger1  the 
next  741215  at  inputs  3  and/or  4,  results  in  it  being  triggered  at  the 
end  of  the  input  pulse.  Thus  its  output  is  delayed  by  the  width  of  the 
input  pulse.  This  fact  was  often  used  in  the  design  of  the  system  now 
described. 

The  pulsing-sampling  rate  of  the  entire  computer  based  PIC  system  is 
control  led,  within  limits,  by  the  SN555  timing  IC.  The  falling  edge  of 
its  positive  output,  is  used  to  trigger  IC-2,  an  SN74121  monostable  multi- 
vibrator. (Note  all  the  remaining  designated  ICs  art  SN74121.)  IC-2 
provides  the  necessary  variable  Lime  delay  between  the  application  of  HV 
and  triggering  of  the  cathode  impedance  rircuitry.  The  trailing  edge  of 
the  positive  output  pulse  of  IC-2,  whose  width  is  controlled  by  P-2, 
serves  to  trigger  IC-3.  The  positive  output  of  IC-3  is  then  fed  into  the 
driving  transistor  Q-l ,  which  inverts  the  pulse  and  drives  the  high 
voltage  transistors  Q-2  through  0-5,  rapidly  to  their  low  impedance  state, 
thus  applying  voltage  to  the  outer  shell  of  the  chamber. 

The  leading  edge  of  the  negative  output  of  IC-2  triggers  IC-4.  Note 
that  this  occurs  sometime  before  the  high  voltage  is  applied.  IC-4 
provides  the  necessary  variable  time  of  occurrence  required  so  that  the 
pulse  which  forces  the  cathode  switching  transistor  into  a  high  impedance 
state  may  be  positioned  to  occur  immediately  after  the  chamber  has  been 
charged  to  its  lull  potential.  The  trailing  edge  of  the  positive  output 
pulse  of  IC-4,  whose  width  is  controlled  by  P-4,  serves  to  trigger  IC-5. 
The  positive  output  of  IC-5  is  then  fed  into  the  driving  transistors  Q-6 
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and  7.  These  invert  the  puise  and  drive  the  two  switching  transistors  Q- 
8  and  9,  forcing  then;  into  a  high  impedance  K20M&)  state. 

On  examining  the  circuit  closely:  one  r.otes  that  the  power  supply 
for  Q~l  and  the  IC's  '!  ,  2,  and  3,  is  floating  shout  350^.  The  reason  for 
this  lies  in  the  way  the  high  voltage  transistors  Q-2  through  5  are  used. 
In  particular,  no  La  that  to  get  Q-2  into  a  r.on-conducting  state  the  base 
must  be  at  a  higher  potential  than  the  emitter,  rfhich  is  at  350  volts. 
To  achieve  this  using  s.  5  volt  pulsing  system,  one  has  simply  to  float 
that  pulsing  sytem  about  350  volts.  Thus  a  voltage  of  352  or-  347  volts 
c&h   be  placed  on  the  base  of  Q-5  forcing  it  off  or  on,  respectively. 

The  critical  parts  to  this  unit  are  the  high  voltage  and  the  cathode 
impedance  switching  transistors.  The  ones  used  in  both  cases  are  state- 
of-the-art  and  were  selected  only  after  an  extensive  search.  The  cath- 
ode impedance  switching  transistors  had  to  have  the  capability  of  con- 
ducting a  large  current  {■-  .5  amp)  at  low  voltage  for  a  short  period  of 
time  ("-  500  nsec)  and  yet  return  to  a  high  impedance  state  in  --  100  nsec. 
Although  this  is  out  of  specification  for  almost  all  transistors  due  a 
basic  design  problem,  the  2N2857  performed  perfectly  vr\6er   these  condi- 
tions. The  high  voltage  transistors  had  to  meet  the  same  basic  require- 
ment as  those  just  listed,  except,  in  addition,  they  had  to  be  able  to 
■isolate  350  volts.  The  2N3743,  although  not  as  fast  as  desired,  performed 
quite  well.  The  KVP  could  drive  a  360  pf  capacitive  chamber-cable  load 
to  350  volts  in  SCO  nsec  and  switch  from  low  to  high  impedance  in  50  nsec. 

The  Raw  Data  Acquisition  system  (RDA)  hdc.   two  basic  tasks;  measure 
the  PJC  output  voltage  v(t  )  and  the  steady  state  current  Igs  for  the  two 
chambers  and  present  the  results  to  the  computer  interface  in  binary 
coded  decimal  (BCD)  format.  The  main  components  of  the  RDA  system  depicted 


in  figure  3-6  are;  2   OL'I  track  and  holes  (T-H),  two  analog  to  digital  (A- 
D)  converters,  a  419  Keithly  picoammeter.  two  318  operational  amplifiers 
(op  amps),  and  four  302  op  amps.  The  logic  of  this  system  is  given 
below, 

The  v(t  )  signal  coining  from  the  cathode  of  the  chamber,  when  the 
system  is  in  the  PIC  mode,  is  Ted  into  a  302  op  amp  with  a  gain  of  one. 
The  signal  is  then  attenuated  by  a  factor  of  2  \x\   gam,  by  using  a  resis- 
tor network,  and  fed  into  the  positive  differential  "input  of  a  318  op 
amp.  The  signa"  is  reduced  by  x   facto*  of  two  because  of  the  fact  that 
the  318  op  amp  is  much  more  stable  when  operated  at  a  gain  of  approxi- 
mately 2.  The  use  of  the  318  op  amp  allows  one  to  apply  a  zeroing  d-c 
shift  to  the  input  as  well  as  to  adjust  the  gain  so  that  the  entire  PIC 
system  has  a  gain  of  one.  The  signal  from  the  318  op  amp  then  serves  as 
the  input  of  the  OEI  5892  track  and  hold  (14!).  The  T-H  takes  a  400  nsec 
sample  at  a  given  time  after  the  collection  voltage  is  applied  to  the 
chambers.  The  time  at  which  the  sample  ^s  taken  is  set  by  the  experimen- 
ter by  adjusting  P-6  of  IC-6.  Note  that  IC-6  is  tripped  by  the  negative 
going  pulse  of  IC-5  and  thus  for  given  pot  settings,  has  a  definite  time 
relationship  with  respect  to  the  time  of  application  for  the  collection 
potential.  IC-7  serves  as  the  T-H  sample  pulse  width  controller. 

The  output  of  the  T-H  is  then  fed  into  the  first  analog  to  digital 
converter  (A~Q1 )  which  presents  the  data  to  the  computer  interface  in  the 
required  BCD  format. 

Once  both  data  channels  have  been  read  in  the  PxC  mode,  the  computer 
then  switches  the  mode  relays  so  that  a  constant  d-c  collection  potential 
is  applied  to  the  chambers  and  the  cathode  of  each  chamber  is,  in  its 
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time,  fed  into  the  419  picoammeter.  The  output,  of  the  ammeter-  is  fed 
v.Vi.0   A-01  which  is  read  by  the  computer,   fhfe  computer  examines  the 
reading  to  determine  whether-  the  ammeter  is  on  the  right  scale.  If  it  is 
not,  the  computer  raises  or  Sowers  the  ammeter  range  by   one  decade  and 
reexamines  the  resulting  419  output.  This  is  done  until  the  meter  reading 
indicates  the  ammeter  has  the  right  range  setting.  This  data  are  then 
stored  and  the  computer  switches  to  the  second  channel  and  repeats  the 
above  sequence  to  obtain  a  correct  current  reading.  Thus  data  are  is  recorded 
in  both  PIC  and  steady  state  modes  automatically.  Once  this  is  done  the 
chambers  are   exposed  to  a  different  intensity  radiation  field  by  moving 
the  chambers  closer  or  farther  from  the  source,  and  the  computer  is  then 
given  a  command  to  repeat  the  above  sequence. 

The  channel  and  mode  switching  system  (CMS)  is  composed  of  transis- 
tors reed  relays,  and  diodes.  It  is  an  integral  part  of  both  the  HVP  and 
RDS  system  and,  as  such,  is  included  in  their  schematics.  Switching  from 
PIC  to  steady  state  mode  is  accomplished  when  the  computer,  through  the 
interface,  provides  a  positive  3  volt  output  which  is  applied  to  the  base 
of  Q--10  and  Q-ll  or  0-12.  This  results  ir:  relays  Rl  and  2   or  3  closing. 
Note  R2  is  closed  when  the  current  from  channel  1  is  being  measured  and  R.3 
is  closed  for  measurement  of  the  current  from  channel  2.  Mote  also  that, 
when  the  current  is  being  measured  in  each  channel,  the  PIC  mode  cathode 
circuitry  is  disconnected.  The  coaxial  switch  serves  to  route  data  to  A- 
Dl  as  described  below. 

R-6  is  closed  in  the  PIC  mode  for  measurment  of  v(t  )  from  chamber 
1.  The  v(t  )  signal  from  chamber  2  is  routed  through  A-D2  when  the  system 
is  in  the  PIC  mode,  R-5  is  closed  as  well  as  R-l  and  R-2  when  measuring 


the  I   from  chamber  1.  R--4  is  closed  as  well  as  R-1  and  R-3  when  measuring 

I  ,  from  chamber  2 , 
ss 

The  timing  sequence  for  one  data  recording  cycle  in  the  PIC  mode  begins 
when  the  cathode  impedance  circuitry  is  tripped  to  its  high  state  by  IC- 
4.  The  positive  going  pulse  of  IC-4  is  fed  into  IC-8.  The  negative 
output  of  IC-8  is  then  inverted  by  Q-21  and  nhen  fed  to  I C - 9  and  IC-10.  The 
pulse  width  of  IC-8  is  set  such  that  it  is  slightly  greater  in  width  than 
that  of  those  used  in  the  T-H  triggering  circuitry  from  IC-6  and  7.  This 
is  to  ensure  that  the  A-D  converters  are  commanded  to  read  their  inputs 
only  after  the  T-H  has  the  data  sample  ready  for  them.  IC-9  and  IC-10 
provide  the  pulses  necessary  to  activate  tne  A-D's.  Since  the  A-D  con- 
verters require  250  msec  to  sample,  the  output  pulses  of  IC-9  and  IC-10 
are  set  accordingly.  l'C-11,  which  is  triggered  by  IC-10  at;  the  end  of 
its  pulse,  indicates  to  the  computer,  through  the  interface,  that  the 
data  is  ready  to  be  read.  The  IC-'M  pulse  width  is  set  so  as  to  allow 
the  HP9S21  2->0  msec  to  record  the  data  from  both  A-D  converters. 

When  the  system  is  in  the  steady  state  mode,  only  A-Dl  is  read  by 
the  computer,  as  described  earlier.  The  codes  used  by  the  computer  to 
control  the  system  and  output  the  resulting  data  are  given  in  the  Appendix. 
The  system  operated  flawlessly  over  the  3  months  of  evaluation.  No 
major  deficiencies  were  found.  The  systems  response  remained  both  linear 
and  free  of  gain  drift  with  the  calibrating  accuracy  of  +2%. 


Detector  Sel ection 
The  fundamental  part  of  any  nuclear  radiation  measuring  system  is 
the  detector.  The  capability  of  a  system,  even  with  the  most  sophisticated 
electronics,  is  ultimately  determined  by  the  sensor.   The  only  detector 
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proven  capable  of  approaching  the  limit?  set  forth  in  the  beginning  of 
this  chapter  is  the   U  fission  chambpr,  It  gives  the  maximum  possible 
gamma  discrimination,  with  relatively  low  burnup  and  has  been  proven  to 
operate  under  the  most  adverse  conditions,  all  of  which  are  necessary 
characteristics  fur  the  PIC  application.  For  these  reasons  two  matched 
pairs  of  RSN-34A-M1  fission  chambers  were  obtained.  A  scale  diagram  is 
shown  in  Figure  3-7.  Of  each  pair,  one  had  a  coating  of  93%  ""U.  Thus 
one  chamber  of  each  set  was  sensitive  to  both  neutrons  and  gammas,  while 
the  other  was  sensitive  only  to  gammas.  As  a  consequence  they  formed  a 
gamma  compensating  pair.  One  set  was  filled  at  the  manufacturing  facility 
with  1  atmosphere  (STP)  of  a  high  purity  Ar-5%N«  gas  mixture.  The  other 
set  was  ordered  with  fill  tubes  attached.  They  were  pumped  down  to  2  x 
10~8  "err  at  300°C,  cooled,  and  then  filled  to  2   atmospheres  with  research 
grade  high  purity  neon.  The  argon-nitrogen  cas  mixture  was  used  because 
of  its  proven  characteristics  in  conventional  fissiun  chamber  operation. 
Neon  on  the  other  hand,  was  reported  to  have  a  volume  recombination 
coefficient  that  vas  independent  of  oernperature;  a  dosired  gas  charac- 
teristic of  the  PIC  reactor  power  measurement  system. 

These  chambers,  along  with  the  previously  described  PIC  system, 
were  then  evaluated.  Chapter  IV  describes  the  experimental  procedures 
used  to  accomplish  this,  as  well  as  the  results  obtained. 
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CHAPTER  W 

EXPERIMENTAL  PROCEDURES,  RESULTS  AND 

METHODS  OF  DATA  PROCFS3TNO 

Initial  S,ystem  Setup 

The  limiting  factor  of  the  sampling  rate  of  the  system  described  in 
the  previous  chapter  was  the  analog  co  digital  conversion  speed  of  com- 
ponents used.  One  complete  data  reading  required  250  msec.  Because  of 
this  limit  the  system's  pulse  rare  was  set  at  1,55/sec. 

The  high  voltage  pulse  width  was  set  at  15  msec,  which,  was  more 

than  sufficient  to  sweep  the  chambers  free  cf  the  ions  '.-'hicri  constituted 

the  steady  state  charge  density.  The  collection  time  for  the  Ne  ions 

in  the  RSN-34A-M1  chamber,  with  a  200  volt:  collection  potential  applied, 

was  100  usee.  For  the  Ar-N^  chambers,  with  a  300V  collection  potential 

i- 

applied,  the  "ion  collection  time  was  found  to  be  130  usee.  A  time  of 
630  msec  was  allowed  for  the  Wiization  in  trie  chamber  to  reach  its 
equilibrium  value.  This  was  found  to  be  more  than  sufficient,  since  the 
ionization  density  reaches  its  steady  state  value  in  less  than  100  msec. 

The  maximum  possible  collection  potentials  were  used  in  each  case 
to  ensure  rapid  collection  cf  the  ions.  The  neon  filled  chambers  had 
200  volts  applied,  because,  above  235  volts,  their  response  was  observed 
to  suffer  from  breakdown.  The  Ar~N?  filled  chambers  were  operated  at 
350  volts,  because  as  stated  earlier,  that  was  the  upper  limit  of  the 
HVP . 

The  system  was  run  continuously  for  two  months,  with  linearity,  gain 
and  zero  drift  being  monitored  periodically  during  that  time.  The  lirearit. 


<";  ■] 
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and  gain  were  checked  over  the  full  voltage  and  current  range  of  interest 
by  "injecting  known  signals.  Over  the  entire  period  of  operation  the  sys- 
tem's response  remained  both  linear  and  at  a  constant  gain  within  a  +2% 
accuracy.  On  the  other  hand,  zero  drift  in  the  PIC  mode  was  found  to  be 
a  strong  function  o\   the  ambient  temperature  of  the  HVP.  The  impedance 
switching  transistors  0-8  and  0-9  were  identified  as  the  cause.  Their 
temperature  coefficient  was  approximately  10  mV/°C.  For  the  measurements 
taken,  baseline  drift  was  carefully  monitored  to  ensure  false  readings 
were  not  obtained. 

Expe r jmen_ta  1  Measurements  &  Results 

Gamma  Compensation 

After  carefully  ascertaining  the  system's  linearity  and  verifying 

its  calibration,  measurements  of  exposure  rate  versus  v(0  were  begun. 

c 

.  60 
The  radiation  source  utilized  was  the  Westinghouse  Han ford  230  kCi   Co 

irradiation  facility,  which  is  depicted  in  figure  4-1.  A  plot  of  the 

radiation  intensity  versus  vertical  position  for  the  port  used  is  shown 

in  figure  4-2,  The  data  points  were  obtained  using  an  RSG-8A  gamma 

ionization  chamber.  The  curve  represents  a  seven-.h  order  polynominal 

fit  of  the  data; 

log(R/hr;   -  6.226   ■■   .4033*  -    .'i939x2  +     02524x3 

-  2.2067xl0"3X4   J-    !.682>:'i0"4Yh  -   3.3963xi0"DAb 

-8  7 
<•  4.0521x10     x   * 

where  x  is  the  digital  position  of  the  chambers  in  the  irradiation  port. 

x  ranges  from  1  to  22   Note  this  curve  fit  equation  was  used  in  the  PIC 

System  Data  Recording  Code  given  in  the  Appendix. 

The  first  system  response  to  be  examined  was  its  ability  to  gamma 

compensate.  To  do  this  the  matched  pair  of  Ar-N?  chambers  were  utilized. 
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Figure  4-1.     Westinghouse  Hanford  230  kCi     "Co  irradiation  facility. 
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They  were  placed  side  by  side  in  a  cylindrical  styrofoam  mold  and  lowered 
in  6  inch  increments  down  into  the  radiation  port,  while  the  computer 
system,  recorded  v(t  )  and  J   of  both  chambers  at   each  point.  The  code 
used  to  control  the  system  is  given  in  the  Appendix  under  the  title  PIC 
System  Data  Recording  Code.  Note  this  code  was  used  to  take  all  the 
v(t  ),  R/hr,  and  I   data  reported  in  this  chapter.  Typical  plots  of  the 
v(t  ),  I  ,  and  R/hr  data  taken  are  shewn  in  figures  4-3,  4-4,  and  4-5. 
It  is  immediately  apparent  from  figures  4-3  and  4-4  that  the  chamber  with 
the   U  deposited  on  its  walls  was  3  times  more  sensitive  to  gamma 
radiation.  This  fact  shows  that  the  composition  of  the  wall  material  has 
a  strong  effect  on  the  gemma  sensitivity  of  a  giver,  chamber.  The  non- 
linearity  observed  in  the  low-end  response  of  the  neutron  sensitive 

234 
chamber  was  dv,e   to  ionization  produced  by  the  °  U  alpha  activity  in  the 

uranium  coating  on  its  inner  walls. 

In  the  plot  of  v(t  )  versus  I   it  can  be  seen  that  the  curves  for 

the  two  chambers  lie  virtually  on  top  of  each  other.   It  shall  be  remembered 

that  S,  the  ionization  source  strength  is  linearly  related  to  Is_,  i.e., 

S  .  13* 

eu 

where  e  and  U  denote,  the  unit  of  electronic  charge  in  coulombs  and 
sensitive  volume  of  the  chamber,  respectively. 

Thus  figure  4-5  shows  that.,  at  a  given  ion  production  rate  S,  regard- 
less of  which  chamber  is  used,  the  resulting  measured  v(tc)  is  the  same. 
This  proves  that  the  chambers  and  electronics  fcr  each  channel  are  closely 
matched,  a  very  basic  requirement  for  gamma  compensation.  Note,  however. 
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as  stated  earlier,  for  a  given  exposure  rote,  the  ionization  rates  in 
both  chambers  are  not  equal,  due  to  the  effect  of  the  uranium  wall  coat- 
ing in  the  neutron  sensitive  chamber. 

Examining  figure  4-4,  one  can  see  the  presence  of  the  linear  and 
second  order  response  regions  of  v(t ,)  versus  R/hr  as  predicted  by  the 
theory  given  earlier.  However,  the  presence  of  transition  regions  and 

their  characteristics  should  be  noted.  In  the  fission  chamber  curve, 

4  5 

this  region,  extends  from  approximately  10  to  at  least  b  x  10  R/hr.  The 

gamma  chamber  appears  to  have  a  much  more  defined  transition  region:  from 

b  6 

10  to  10  R/hr.   fhese  effects  are  observed  in  figure  4-5  as  well.  The 

broad  exposure  range  covered  by  the  transition  region,  coupled  with  the 

fact  that  this  region  of  response  does  not  coincide  between  chambers, 

makes  gamma  compensation  difficult.  The  gamma  chamber  signal  cannot  be 

simply  multiplied  by  a  constant  and  subtracted  from  the  fission  chamber 

to  give  full  compensation  over  the  full  range  of  exposure  rates,  because, 

as  seen,  due  to  the  increased  gamma  sensitivity  of  the  fission  chamber, 

the  two  chamber  responses  are  not  linearly  related..  Thus,  to  accomplish 

gamma  compensation,  using  the  chamber  selected,  more  than  simple  differ- 

7 
ential  gain  controlled  inputs,  as  were  used  by  Cooper  and  suggested  by 

9 
Ellis,  were  required.  To  compensate  using  nonlinear  electronics  would 

have  been  complicated.  On  the  ether  hand,  compensating  through  computer 
methods  appeared  relatively  simple  and  straight  forward. 

Computer  based  compensation  was  accomplished  as  follows.  The  cham- 
ber responses  shown  in  figure  4-4  were  curve  fitted  to  eighth  order  poly- 
nomials, using  Chebyshev  Pclynominals.    A  description  of  the  code  is 
contained  in  the  Appendix.  The  resulting  equations  were; 


X  =  B.I 6!  +  i.4986y   .4590y2  -  .2C15y3  -  .4586y4  (4-1) 

+  ,2038y5  ■*  .366:;/  +  .0S162y7  -  .009803y8, 
X  =   5.6355  +  1 .532: z  +  .'I500z2  -  .3241z3  +  .0488424         (4-2) 
+  .41387°  +  .251 9z6  +  .05276::7  +  .0025762", 
where;   x  "  1  on  R/hr, 
y  =  log  vf(tc), 

z  -"■  log  v  (tc), 

Vjr(t  )  -  the  PiC  voltage  signal  from  the  fission  chamber, 

v  (t  )  -  the  PIC  voltage  signal  from  the  gamma  chamber. 

A  code  was  then  written  for  the  HP9821  which  caused  it  to  record  the  v(t  ) 

c 

signals  from  the  chambers,  at  given  exposure  rates,  compute  the  measured 
R/hr  from  equation  4-1  or  4-2  for  the  corresponding  chamber,  subtract  the 
results,  compute  the  compensation  error  and  finally  output  all  results. 
The  code  is  listed  in  the  Appendix  under  the  title  of  Compensation  Code. 

Table  4-1  contains  a  set  of  typical  results  of  this  application. 
The  sixth  column  in  the  table  plotted  in  figure  4-6  indicates  the  com- 
pensation results.  The  first  siy  values  are  as  large  as  they  are  for  two 
reasons.  The  first  is  that  the  curve   fit  did  not  fit  the  low  data  points 
well  and,  second,  the  system's  electronic  stability  was  +2mV,  which,  as  one 
can  see,  has  a  large  effect  on  the  computed  R/hr  values  at  the  low  ex- 
posure rates.  The  remainder  of  the  values  in  that  column,  nevertheless, 
indicate  that,  at  least  for  fixed  temperatures,  reasonable  compensation 
can  be  obtained.  The  fluxuations  that  do  exist  are  due  to  curve  fitting, 
chamber  positioning,  and  the  precision  of  the  measuring  system.  The  last 
of  these  was  measured  at  +1%.  The  other  two  were  difficult  to  accurately 
determine,  but  combined,  are  on  the  order  of  +5%. 
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Figure  4-6.  Gamma  compensation  as  a  function  of  exposure  rate. 


PIC  Chamb er  Temperat \rr e  Respon s e 

Once  gamma  compensation  was  proven  feasible,  using  the  computer 
based  method,  v(t  )  for  both  chamber  sets  was  examined  as  a  function  of 
temperature  at  various  exposure  rates.  This  was  accomplished  by  con- 
structing the  3.8  inch  diameter,  double  tubed  oven  shown  in  figure  4-7. 
The  Ar-Np  chambers  were  inserted  into  the  oven  and  connected  to  the  long 
RG62-U  leads  from  the  PIC  HVP  by  2-foot  lengths  of  high  temperature 
mineral  (SiOp)  insulated  cables.  However,  when  the  Ne-filled  chamber 
measurements  were  taken,  air  coax  was  used  in  place  of  the  mineral  cables 
because  they  suffered  from  voltage  breakdown  at  400°C.  The  oven  was 
capable  of  reaching  temperatures  up  to  550°C  when  used  in  conjunction 
with  a  120  volt  variac. 

The  cylindrical  oven-chamber  assembly  was  placed  in  the  same  port  as 
was  used  for  the  gamma  compensation  experiment.  The  data  were  then  taken 
for  both  the  Ar-N„  and  Ne  chamber  pairs.  The  most  pertinent  of  this  is 
displayed  in  figures  4-8  through  4-23. 

Some  basic  chamber  response  characteristics  are  immediately  apparent. 
The  most  obvious  being  the  rise  in  leakage  current  as  the  temperature  is 
increased.  This  fact  accounts  for  the  drastic  change  in  the  slope  of  the 
beginning  data  points  as  the  temperature  is  increased.  Note  that  the  I 
values  became  next  to  useless  as  the  temperature  was  increased,  because 
the  leakage  current,  far  exceeds  the  current  caused  by  the  radiation.  On 
the  other  hand  v(t  },  even  at  the  highest  temperature-,  maintains  its 
basic  functional  response  tc  the  exposure  rate.  {Note,  in  figures  4-7 
through  4-22  the  tabulated  values  from  top  to  bottom  correspond  to  ver- 
tical port  positions  14,  13,  12,  11.5,  11,  10.5.  10,  9.5,  9,  8.5,  8,  7.5, 
7,  6.5,  6,  5.5,  5,  4.5,  4,  3-5,  3,  2.5,  end  1.5,  in  that  order.) 
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Figure  4-7.  Cross  sectional  vie\>   of  the  over;  b\ 
matched  pair  of  RSN-34-M.1  chambers 
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To  glean  information  from  the  data  more  easily.,  the  leakage  current 
was  subtracted  from  all  current  values  arid  the  ratio  of  I  /v(t  )  was 

SS   '  c 

computed  for  all  data  point:;.  Tables  4-2  and  4-3  contain  the  results 

for  the  Ar-Np  and  Ne  chambers,  respectively.  The  blank  spaces  indicate 
the  result  obtained  was  meaningless.  Figures  4-24  and  4-25  present  the 
tabulated  Iss/v(t  )  ratio  at  position  3  as  a  function  of  temperature  for 
He  and  Ar-N?  data,  respectively. 

The  reason  for  calculating  the  ratio  was  to  smooth  out  fluctuations 
in  the  data  caused  by  vertical  and  rotational  positioning  differences 
between  data  sets.  Since  I   is  constant  for  a  given  exposure  rate, 
over  the  temperature  range  of  interest,  v(t  )  accounts  for  the  variation 
in  the  ratio  *ss/v(tc).  If  v(tc)  gets  larger  in  the  second  order  region, 
this  indicates  that  the  volume  recombination  coefficient  is  decreasing. 
Note,  as  we'll,  that  as  v(t  )  nets  larger  I.  /v(t  )  gets  smaller. 

v..  iS     c 

With  these  facts  in  mind  examine  the  Ne  data  in  Table  4-2  and 
Figure  4-24.  In  particular  choose  a  position  where  v(t  )  is  in  the 
second  order  region,  for  instance,  position  3.  The  ratio  of  I  /v(t  ) 
changes  considerably  over  the  temperature  range  of  interest.  At  25°C 
the  ratio  for  the  gamma  chamber  is  4.36  x  10"°  while  at  477°C  it  is  1.32 
x  10  .  The  same  temperature  range  for  the  fission  chamber  gives  a 
ratio  range  of  4.34  x  10"'°  to  1.87  x  10  .  Note  that  the  relative 
changes  are  approximately  equal  for  both  chambers.  As  one  will  see  this 
is  not  true  for  the  Ar-N?  chambers.  This  indicates  the  a?   decreases 
with  increasing  temperature  for  Ne;  a  fact  which  appears  to  contradict 
Sanders  results.  Sanders  found  u,.  constant  for  No  from  25°C  to  300°C. 

L 

Examining  the  data  given  in  the  lower  part  of  Table  4-2  (positions  1.5, 
2.5,  3,  3.5)  more  closely,  one  notes  that  the  ratio  is  virtually  constant 
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up  to  232°C.  At  28?.°C  it  takes  a  sudden  drop  a^d   continues  down  over 
the  rest  of  the  temperature  range.  Thus  rather  then  contradicting 
Sanders  data,  it  validates  his  conclusion  to  some  extent  that  a?  for  Ne 
is  essentially  constant  from  25°r  to  approximately  300°C.  Note  also 
that  a  1  atmosphere  chamber  was  used  in  Sanders  experiment  while  the 
ones  used  here  contained  Z   atmospheres  of  Ne. 

The  gradual  decrease  in  the  ratio  at  a  given  temperature  is  due  to 
the  transition  between  weak  linear  diffusional  losses  and  stronger 
volume  recombi national  losses  as  the  steady  state  ion  density  increases. 

Table  4-3  contains  the  Ar-N2  data.  For  the  highest  exposure  rates, 
where  volume  recombination  predominates,  the  l  /v{t  )  ratio  changes  by  a 
factor  of  10  for  the  gamma  chamber  and  by  a  factor  of  5  for  the  fission 
chamber.  There  is  no  theoretical  explanation  for  this  if  both  chambers 
were  prepared  in  the  same  manner  and  filled  with  the  same  gas.  A  possi- 
ble explanation  is  that  these  chambers  were  not  thoroughly  baked  and 
pumped  before  filling.  The  fission  chamber,  having  its  walls  coated  with 
pure  uranium,  maintained  a  purer  gas  mixture  as  the  temperature  was 
raised  because  the  gas  impurities  in  the  walls  were  trapped  by  the  coat- 
ing. The  gamma  chamber,  however,  having  no  such  protection,  had  wall 
trapped  gas  impurities  mixed  in  with  the  Ar-M,.  fill  gas.  Such  impurities 
would,  in  general,  enhance  the  temperature  dependence  of  a**   an  effect 
which  is  apparent  in  the  data. 

Little  was  said  of  the  transition  and  first  order  data  regions  for 
two  reasons.  The  first  being,  it  is  a  well  known  fact  that  diffusion, 
which  controls  losses  in  the  first  order  region,  is  both  temperature  and 
pressure  dependent  in  rather  complex  ways.lf   thus  a  simple  explanation 
of  the  results  does  not  exist  and  it  would  be  inappropriate  to  go  into  it 


in  this  paper,  since  it  is  not  needed  and  man}  excellent  works  on  diffusion 
have  been  published.  The  second  reason  is  that  for  practical  reactor 
operation  gamma  compensation  is  not  necessary  until  an  exposure  rate  of 
greater  than  1(T  R/hr  is  experienced.  As  one  can  see  in  figure  4-4, 

r 

above  10  R/hr  second  order  effects  predominate. 
Sou rce s  of  Errors 

The  main  source  of  error  in  this  experiment  was  detector  position- 
ing. A  +5%   error  in  positioning  produced  as  much  as  a  +]5%   error  in  the 
measured  signal.  The  reason  for  the  magnification  was  the  rapid  vari- 
ation of  flux  with  position.  In  the  temperature  measurements  the  posi- 
tional error  became  even  more  aggravated  due  to  the  oven.  The  oven  was 
allowed  some  rotational  freedom  as  it  was  moved  vertically.  Thus,  the 
chambers  were  rotated  in  the  horizontal  plane  with  respect  to  the  source. 
This  effect  was  minimized  when  the  styrofoam  mold  was  used,  but  when  the 
oven  was  used,  appreciable  (+_5%)  exposure  rate  fluctuations  were  experi- 
enced at  the  chambers  due  to  self  shielding.  This  variation  accounts,  in 
part,  for  the  apparent  staggered  look  of  the  data  tabulated  in  Table 
4-2  and  4-3.  Only  trends  are  meaningful  in  this  data  and  point  by  point 
comparison  should  be  avoided. 

Another  source  of  error  in  measuring  v(t  )  was  the  drift  in  the  base 
line  caused  by  the  temperature  dependence  of  the  impedance  switching 
transistors  (10  mV/°C).  This  error  was  almost  impossible  to  accurately 
determine.  However,  the  base  line  was  monitored  before  and  after  each 
set  of  measurements,  if  it  was  found  to  have  drifted  more  than  +  10  rnV 
the  set  of  data  was  discarded  and  the  measurement  repeated. 

As  stated  in  the  first  pav"t  of  the  chapter  the  electronics  remained 
both  linear  and  accurate  to  within  +2%,   over  the  3  months  the  system 
was  used. 
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In  measuring  v(t  )  a  source  of  error  other  than  instrument  error  was 
present,  I  ,  which  contributes  to  v(t_)  only  minutely  at:  low  exposure 

S  5  C 

rates,  at  the  high  exposure  rates  could  begin  to  significantly  contribute 
to  the  v(t  )  measured.  Ideally  v(t  )  should  ba  constituted  of  only  the 
steady  state  ion  density  collected.  The  equation  which  describes  the  I 
current  contribution  is, 

vc..(t)  =  I  R  (I  -  exp(-t/RC)), 

where  RC  is  the  circuit  time  constant.  For  t  «   RC  this  reduces  to 

Note  that  v„.(t)  can  be  reduced  by  reducing  R,  just  as  long  as  t  <<  RC. 
b  s 

Also,  to  avoid  affecting  v(t  )  adversely  RC  must  remain  significantly 

greater  than  t  ,  In  this  experiment  R  was  20  Mft,  C  was  £00  pf  and  t  was 
3  c  r      c 

IOC  usee.  Thus  RC  was  8  msec,  a  factor  of  3000  greater  than  t  .  At  the 
high  exposure  rates  v  (100  usee)  constituted  approximately  5%  of  the 
measured  v(t  ).  This  could  have  been  reduced  by  a  factor  of  approximately 
20  by  reducing  R  to  1  Ha.     Such  a  reduction  would  have  had  a  minimal 
effect  on  the  true  value  v(t  ). 

Another  possible  source  of  error  incurred  during  the  neon  chamber 
temperature  measurements  involved  the  2  foot  air  coax  used  to  connect  the 
chamber  to  the  RG62U  cables  coming  from  the  HVP.  Since  approximately  a 
two  inch  section  of  the  air  coax  was  partially  heated  during  the  experi- 
ment the  chance  in  the  characteristics  of  the  air  coax  with  temperature 
could  possibly  have  had  a  measurable  effect  on  the  apparent  response  of 
the  neon  chamber.  However,  it  is  felt  that  6ue   to  the  gas  kinetic  charac- 
teristics of  air  and  the  relatively  small  volume  involved,  the  air  coax 
temperature  effect  on  the  measured  signal  was  negligible. 


CHAPTER  V 
CONCLUSIONS 


To  verify  the  Pulsed  Ion  Chamber  applicability  to  the  measuring  of 
the  full  range  of  reactor  power,  three  basic  steps  were  required,  A 
computer  based,  solid  state,  dual  chamber  PIC  system  had  to  be  designed 
and  built.  Gamma  compensation  for  the  expected  mixed  order  ("linear 
and/or  second  order)  response  region  had  to  he  clearly  verified.  The  PIC 
mode  chamber  response  had  to  be  proven  to  dp  independent  of  temperature 
over  a  range  of  25  to  500°C. 

The  design  and  building  of  the  computer  based  PIC  system  took  con- 
siderable time  and  effort,  but  the  end  result  was  success.  It  met  all 
the  design  requirements  and  more.  Its  versatility  simply  as  a  plasma 
diagnostic  tool  is  apparent.  One  has  full  control  over  all  pulsing- 
sampling  sequences,  as  one  can  see  in  figure  3-3  all  the  timing  sequence 
potentiometers  are  on  the  front  panel.  The  system  is  free  from  the  time 
jitter  and  position  orientation  problems  which  plagued  the  most  advanced 
mercury  wetted  reed  pulsers  previously  used.  The  impedance  switching 
circuit  has  a  response  time  on  the  order  of  50  nsec,  due  to  the  use  of 
the  most  advanced  gigahertz  switching  transistors.  This  in  itself  is  an 
improvement  of  a  factor  of  5  over  the  previous  systems.  Thus  the  system 
is  a  significant  advancement  of  the  PIC  state  of  the  art. 

Verifying  the  feasibility  of  gamma  compensation  was  accomplished 
relatively  easily  using  the  PIC  computer  based  system.  Compensation,  at 


<-e3 


a  given  temperature  for  the  exposure  rates  where  compensation  would  be 
required,  can  be  obtained  to  within  +22  or  the  total  gamma  signal, 

The  FIC  mode  chamber  response  was  proven  to  be  temperature  dependenl 
over  the  range  of  25  to  5G0CC  for  Ne  as  well  as  the  Ar-N?  fill  gas 
mixture,  however,  as  stated  earlier,  the  Ne  response  was  constant  up  to 
'\.282°C!  this  constant  range  could  be  extended  by  reducing  the  pressure 
of  the  Ne  fill  gas,  The  300°C  temperature  range  is  the  operating  range 
of  current  light  water  reactors.  Thus  this  PIC  system  offers  all  the 
required  characteristics  and  more  for  an  in-core  wide  range  neutron 
measuring  system  for  present  light  water  reactors, 

A  fill  gas  whose  PIC  response  is  independent  of  temperature  over 
the  25°C  to  500°C  range  must  he  found  before  the  system  could  realis- 
tically be  applied  to  power  measurements  in  the  now,  more  advanced  KTGR 
and  LMFBR  reactors.  Seeing  that  this  is  apparently  the  only  remaining 
stuml  lira  block  to  the  use  of  the  PIC  system  for  wide  range  reactor' 
power  measurements  in  these  reactors  a  concentrated  effort  to  find  a 
suitable  gas  should  be  launched. 

In  conclusion,  since  new  findings  and  advancement  or   the  state-of- 
the-art,  both  of  which  have  been  accomplished  here,  are  the  essence  of 
research,  this  was  indeed  a  successful  endeavor. 


APPENDIX 
WP   9821  Codes 

The  following  HP9821  programs  were  'ised  to  control  the  PIC  system, 
redd   and  record  data,  perform  data  analysis,  and  finally  output  the  data 
in  plot  form. 

PIC  System  Data  Recording _Cod_e 

This  code  records  the  I  ,  v(t  )  and  R/hr  data  used  in  this  thesis. 
j  s     c 

0;  PRT  "THIS  PROGRAM" 

1:  PRT  "RECORDS  S,S.;| 

2:     PRT  "  CURRENT  VS.  PIC" 

3.  PRT  "VOLTAGE  SIGNAL" 

PRT  "DATA  AND  STORES" 

PRT  "IT  IN  FILE" 

PRT  "GIVEN," 

SPC  3 

ENT  "WTB  NO.  C--M0DE  1",  R12 

ENT  "WTB  NO.  C-MODE  ?.",  R13 

ENT  "VOLT  PULL  SCALE",  R7 

ENT  "VOLTS  ZERO  SCALE",  R6 

ENT  "I  FULL  SCALE?",  RIO 

ENT  "I  ZERO  SCALE?" ,  R9 
21  ->   RIG;  0  ■>  R17 


8.-, 


8b 


In  the  above  statements  the  expected  value  limits  cf  v(t  )  and  I 

c      ss 

are  input  as  well  as  the  range  control  numbers  for  the  Keith! ey  419  pico- 

amp;  statements  S  and  9. 

15:  WTB  1,  192 

16:  CMQ  "1\I5":   FMT  *;  RED  13,  A 

17:  1  -y  X 

18:  RED  3,  Y,  Z 

19:  X  +  1  ->  X;  IF  x  <  200;  JMP  -1 

20:  CMD  "?V5";  FMT  *;  RED  13,  B 

21:  B-A  ■>  C;  IF  C  >  100;  INT  (C/100)  *  60  +  (C-INT  (C/100)  *  100)  -  C 

22:  C  *  1E6/200  ■*   R5;  FLT  5;  PRT  "PERIOD  IN  10-6  SECONDS  — \  R5 

Statement  15  contains  the  command  which  switches  the  system  into  the 
PIC  mode.  The  remainder  of  this  section  measures  the  period  of  the 
system  utilizing  the  HP  clock  which  is  addresser!  in  statements  15  and  20. 


ENT  "POSITION  ?",  X 

"2I!;  WTC  1,  R12;  R3  +  X 

6.226376  +  .403346  X  -  .1 93881  x  t  4  ■*  Y 

Y  +  1.168232E  -  4X  ■*■   5  -  3.3Q6338E  -  6X  t  6  f  4.052068E  -  8X  +  7  ->  Y 

10  *  Y  -»  R  (RI6  +  4) 

The  vertical  port  position  is  entered  by  the  experimenter  and  the 
R/hr  is  computed  from  the  seventh  order  curve  fit  descirbed  in  Chapter 
IV. 

Statement  24  switches  the  system  into  the  steady  state  mode  in  order 
to  measure  I  ,  for  channel  1(1  ,)• 

S  S  S  S  : 

28;   "3";  GS3  "DELAY" 

29:  RED  3,  X,  Y 

30:   IF  Y  <  200;  IF  R12  <    170;  R12  +  1  -*  R1Z;  GET  "2" 
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31; 
32 

33- 
34: 
35: 
36. 
37: 


40: 
41: 
42: 

43: 


the 

44 
45 
46: 

47: 
48: 
49: 

50' 


IF  Y  >  3000;  IF  R12  >  160;  Ri2  -  1-  ■*  R12  GTO 
GSB  "DELAY" 
RED  3,  X,  Z 


ABS  (Z-Y)  >  60;  GTl 


0 


RED  3,  X,  Y 

Z  +  Y  +  Z;  C  +  1  '■*   C;  IF  C  <  8;  JMP  -  1 
Z/10  ->  1 

1*   .0033333  *  10  t  (158  -  R12)  •*  RR<6 
IF  R12  =  165;  RR1 6  -  .1E-7  -  RR16;  GTO  "4" 
IF  R12  =  166;  RR16  -  .05E  -  8  -  RR16;  GTO  "4" 
IF  R12  =  167;  RR16  *  1E9  ■*■  X 

IF  R12  -  167;  .03948  +  .8529X  +  79G8E  -  6  *  X  +  2  -  4.202E  -  4  *  X 
3-*Y;Y*lE-9  +  RR16 

This  section  measures  I  ,   10  times  and  then  computes  and  stores 
average. 
"4";  WTB  1,  R13 
GSB  "DELAY" 
RED  3,  X,  Y 

IF  Y  <  200;  IF  R13  <_  106;  Rl 3  +  1  ->   R13;  310  "4" 
IF  Y  >  3000;  IF  R13  >  96;  R13  -  1  -►  R13:  GTO  "4" 
GSB  "DELAY: 
RED  3,  X,  I 
IF  ABS  (Y  -  z)  >  60;  GTO  "4" 


53:  RED  3,  X,  Y 

54;  z  +  Y  >  Z;  C  +  1  -*   0;  IF  C  <-  8;  JMP  -  1 


56: 
57: 

58: 

59, 
50: 
5] 
62: 

63: 
64; 
65: 

66: 
67: 
68: 
69; 


Note 

70- 

/l  ; 

72; 
/  5, 

74: 
75 
76 
77- 


78: 


2/10  -*  Z 

Z  *  .003333  *  10   '•   (34    -  R13)    ►  R  (R16  +  2) 

R   (R16  +  2)   ■»-  X 

IF  R13  -  101;   x  ■•    .  1E-7  -  R    (RIG  +  ?};   CTO   "B" 

l!;  R13  =  102;   X  -.05E-8  -   R   (R16  +  2);   GTO   "B" 

IF  R13  1   103:.   GTO   "B" 

X   *  1E9  -*■  X;    .03948  +    .8529X  +  75G8E-C   *  X  +  2  -  4.202E-4  *  X  t  3  ->  Y 

Y  *  1E-9  ■*  R  (R16  +  2) 

This  section  measures  and  store::  the  I   value  for  channel  2. 

"B";  WTB  1,  192 

0  -  C;  0  ■*■  X;  0  ->  Z;  0  ->  Rl  8;  0  ->  R19 

35  -  X;  GSB  "DELAY" 

RED  3,  Y,  X;  RED  5,  Y5  Z 

X  +  R18  ■+  Rl S;  1  +   Rl 9  ->  R19;  C  +  1  ->   C;  IF  C  <  9;  JMP  -1 

RIB/10000  -»  R  (R16  +  1 } ;  R19/1000  -»  R  (R16  +  3) 

FLT  4;  PRT  R  (RI6  +  4),  R  (R16  +1 ) ,  RR16,  R  (R16  +3),  R  (R16  +2}; 

R16  •:-  6  -*•  R16;  SPC  1 

v(t„ }  tor  channels  1  and  2  are  measured  and  stored  in  this  section, 
c 

10  values  of  each  are  treasured  and  then  averaged  together. 

ENT  'MORE?  YES  +  1  NO  -*  0",  Y 

ENT  ''POSITION?11,  R3 

IF  Y  >  0;  GTO  "2" 

ENT  "NE^  SET?  Y  +   1  NO  ■+  0"s  Y 

IF  Y  >  0;  R;6  ■+  RR17;  Rl 7  +  1  ■*  Rl/;  GTO  "2" 

R16  ■*  RR17 

FNT  "TAPE?",  A,"  FILE?',  A 

RCF  A,  R  (R16  -  1) 
GTO  "END'1 


This  section  asks  if  more  data  are  to  be  taken.  If  so,  the  cycle  is 
repeated.  If  not  the  data  dre   stored  on  tape  in  the  specified  file. 
79:   :'DtLAY" 

80:  CMLi  "775";  FM1  *;  RED  13,  A 
81:  "0";  CMD  "?V5";  RED  13,  3 
82:   IF  X  }   35;  10  ■*  X 
83:   IF  B-A  <  X;  GTO  "0" 
84:  RET 
85:  ,:END";  END 

This  is  a  subroutine  addressed  after  switching  from  channel  to 
channel  cr  mode  to  mode.  It  delays  the  data  recording  for  10  to  35  sec 
to  allow  the  system  to  stabilize. 

Curve  Fitt ing  By  Chebys h ev  Polyn omina Is 

This  program   fits  a  least-squares  curve  to  a  set  of  given  data 

points  (x,,  y-j),  (x?,  y? ) (x  ,  y  ,)  where  the  x.  lie  in  an  interval 

(a,  b)  and  are  equally  spaced.  The  user  specifies  a  degree  m  and  the 

program  outputs  the  coefficients  a,,?  a,,...,  a  of  a  polynomial  P(x)  =  a^ 

+  a-,x  +...+a  x  passing  near  or  through  each  incut  point. 
1      m   •  ■ 

The  program  determines  the  a.  by  considering  P(x)  as  a  linear  combin- 
ation of  Chebyshev  polynomials  T.(x),  P(x)  -  c0TQ(x)  +  c-,T.,(x)  +...+ 

c  T  (x)  and  applyinq  the  least-squares  criterion  to  the  expression 

m  m       •  '  -   •'  ' 


as 
v\  -  P(x)    to  give  a  system  of  simultaneous  equations  —  •r  -   0, 


i=l 


cCJ 


j  =  0,  1,...,  m  from  which  the  c.   can  be  determined.  Calculation  of  the 
c.  is  facilitated  by  using  the  orthogonality  properties  of  the  Chebyshev 


90 
)Olynomials  and  evaluating  P{x)  at  special  points 


y  s  £2|4-^L.L  y.l.    within  the  interval  (-1,  1)  to  force  off-diaqonal  ele- 


merits  to  be  zero  in  this  system  of  equations.  Corresponding  values  y. 

are  needed   in  the  system  of  equations  in  order  to  be  able  to  solve  for 

the  c.  The  program  obtains  these  by  applying  a  linear  transformation  to 

the  x.  to  bring  the  x-  within  the  interval  (1,  1)  arid  then  using  these 

transformed  values  and  applying  the  Lagrange  interpolation  formula  to  the 

>;.  to  obtain  the  y..  The  system  of  equations  is  then  easily  solved  for 

the  c_. .  The  program  then  applies  a  linear  transformation  lo  x  to  change 

the  expression  P(x)  =  cQT0(x)  +  cft^X)  +...+  ct[Tr.(x;  to  the  form  P(x)  - 

a„  +  a-,x  +...+a  x™  over  the  original  interval  (a,  b). 
0    I       m 

0:  TBL  4;  ENT  "DEGREE?",  7.  +  1  -*  R6,  "XI?",  R8  -*  X,  "DELTA  X?n,  R9; 

CFG  13;  TBL  2 
i:  PRT  "DEGREE";  SPC  1:  PRT  Z;  SPC  2;  PRT  "X" ,  "Y";  SPC  1;  12  ■*  RO  ->  Z 
2:  ENT  "Y?\  RZ;  IF  FLG  13  -  0;  PR"!"  x,  RZ;  bPC  1 ;  X  +  R9  *  X;  Z  +  1 

-  Z  ->  Rl :  JMP  0 

SPC  1;  PRT  "COEFFICIENTS",  SPC  i 

(((Rl-12  -  R7)  +  R  (l  ->  A)  ->  R2)  +  R6  ■>  R3)  +  R6  ■>  R4 

COS  ((tt/2)  (2A-U/R6)  ■+  R  (R6-A  +  R2);  JMP  (A  +  i  ->■  A)  >  R6 

R7  -  1  ->  A;  2/ A  ->  ■  B;  -1  ->  RR1 ;  Rl  ->-  C 

RC  +  B  *   R  (C-+  1);  JMP  (C  +  1  .->  C)  =  R2  -1 


9; 


14: 

lb: 

)  6 

17 

18 

19: 


23: 

24: 


27: 

2S: 

29: 

30: 

31 

32 

33 

34' 

35: 

36: 

37 


IF  R   (R2  +    \}    >  R   (Rl   +  b).;  B   i-  1   ->■  B;   JMP  0 

B  +  (B  =  1)   -  (B  =  R7  -1)  ->  C 

(R7-I)    (R  (R2  +■  A)   -  R  (Rl   +  C  -l))/2  +  Y 

-Y   (Y  -   1)    (Y  -  2)   R   (R0  +  C   -2  ->  C)/6   +   (YY  -1 )    (Y  -  2)   R   (C  +  l)/2   •>   Z 

Z  -   (T  +  1)    (Y  -  2)   YR   (C  +  2)/2  +  Y   (YY  -    i)   R   (C  +  3)/6  +  R   (R4  +  a) 

IF   R6  ■-:    (A  •!-  1   •>  A);   GTO   -5 


0  ■+  Y;  0  -  B;   IF  A  =  1;  GTO  25 
TF  A  >  1;   GTO  27 

Y  +  R   (R4  +  B)  ■*  Y;   JMP   (8  ■♦    1   ->  B)    >  R5   -1 
GTO  31 

Y  +  R   (R2  +  B)   R   (R4  +  6)  -►  Y;   JMP   (B  +  1  -*■  B)    >  R6  +   1 
GTO  31 

1  +  RR1 

R   (R2  +  B)  ->   R   (Rl   +  1);   2  -*■  C 

2R   (R2  +  B}   R   (Rl    +  C-l )    -R   (Rl    +  C-2)  ■>  R   (Rl   4  C);   JMP   (C  +  1  ->  C)   >  A 

Y  r  R   (R4  +  B)   R   (Rl   -f   A)   -  Y;    IF   (B  +  1   ■*■  B)    <  R6    -1;   GTO  28 
2Y/Ry    y  R    (R3  +  A);    IF   (A  +  1   ->  A)   <  R6   - 1 ;   GTO  21 

RR3/2  ->  RR3 

(R8  •>  C)   4   R9   (R7   -  1)  f  B 

(B  +  C)/2  -  RIO;    (B  -  C)/2  •>  Rll 

RR3  ->■  RR4;   R   (R3   «■  1}   ■+  R   (R4  +   1);   0  -*  C 

IF  R6   <  2;   GTO  46 

0   -   R   (Rl    +  C)   -  R   (RO  +  C)  ->  R   (R4    +  C  +  2);   JMP   (C  +  1   ->  C)    >   R6   -I 


!;2 


3d:     1  -»■  R  (Rl   +  1);  2  ->  C 

39:     C0S-(ttC/2)  -*  RRO;   1   -*  B- 

40:      2R   (Rl   +  B  ~   1)    -R   (RO   f  5}  ->  R   (RO  +  B);   JHP   (B  +  1  ->  B)   >  C 

41 :     0  ->  B 

42:      R   (R4  +  B)  +  R   (R3  +  C)   R   (RO  +  B)   ■>  P   (R4  +  G);   R   (Rl   +  B)   ->  A; 

R  (RO  +  B)  ->  R  (Rl   +  B);  A  -  R   (RO  +  B) 
43,      IF   (B  +   I    ■>  B)  £  C:   GTO  42 

IF  (C  +  1  -*•  C)   <  R6  -  1;  GTO  39 


44 

45 
46 

47 
^3 
49 
50 
5 1 


RR4  ->  RRO;    1   +  C 

0  ->  R  (RO  •)   C);  JMP   (C  +  1  ->  C)  =  R6 

1  ••>  C 

R  (RO  +  C)  +  R  (R4  +  C)/R11    +  C  ->  R  .(RO  +  C) ;  1  -*  A  ->  R5  •>  B;  C  +  1  ->  Y 
B  +  X;   XR5  ->■  R5;   R5R11    !■  C  ->  1\  A   (Y  -  X)  ~  A 

R   (RO  +  C  ••  B)  +  A*R10  +  B*(-l)   I  B*R   (R4  +  C)/Z  ->  R  (RO  +  C  -B);   IF 
(B  <-  1  +  B)   <  C;  GTO  50 
IF  (C  +  1   -»■  C)  f  R6;  GTO  49 
0  ->  A 
PRr  A,   R   (RO  +  A);   SPG;  JHP   (A  +  1   ->  A)   =   R5 

RO  ■*  R7:    PR!    "— ";   SPC  13 

GTO  0 
END 

Conipe nsation  Code 
This  program  was  written  to  determine  iiow  accurately  gamma  compensa- 
tion could  be  achieved  utilizing  digital  computer  methods  rather-  than 
analog  methods . 


0 

:  PRE  "---  DATA  ---";  SPC  ?,  PR!  "R/HR";  SPC  1;  PRT  "V  "3"",  "V  "4"" 

',  SPC  1 

] 

:  PRT  "R  =  F  (V3)'!,  "R  -  F  (V4) ' ;  SPC  1 

£. 

:  PRT  •'%   D1F  of  V3  to  V4";  SPC  1 

> 

;  PRT  "%/DIF.  to  R/HR";  SPC  3 

4 

:  21  -y   R16;  0  ->  R17 

5 

.  WTB  1 ,  192 

6 

"A";  ENT  "POSITION  ?",  X'3  IF  X  <  0;  GTO  "L" 

7 

6.226376  +  .603346X  -  ,193881X1?  +  .02637X+3+Y 

8 

Y  -  .00220673X+4  +  1  .168232E-4X+5  -  3.396338E  -6X+6  +  4.052068E-8X+7+Y 

9 

10  +  Y  -:-  RR16 

10 

0*C;0  +  XjO  ■*   Z;  0  ->   RI  8;  0  ■>  R19 

In  the  above  section  the  experimenter  enters  the  vertical  port  posi- 

t i on  and  the  exposure  rate  is  calculated  and  stored. 

11 

35  ■*  X;  GSB  "DELAY" 

1? 

RED  3,  Y,  X;  RED  5,  Yt  7. 

13 

X  +  R18  -  R18;  Z  +  R19  -*  R19;  C  +  1  ->  Ct  IF  C  <_  9;  JMP  -1 

14 

RI 8/1 000  ->  R  (R16  +  1);  Rl 9/ 10000  -*  R  (R16  ■<  2) 

15 

LOG  R  (R16  +  1)  +  X 

16 

5.1661489  -^   1.4982995X  f  .4589825Xr2  -  .20!5837Xi3  -  .458746X+4  ■*-  Y 

17 

Y  +  .2037761X+5  *  .36659235X^6  +  .08162451 X+7  ■*  Y 

18 

-  9.80343538E-3Xta  +  Y  ->  Y;  lOtY  +  R  (PIS  +  3) 

19 

LOG  R  ( P.I  6  +  2)  ■+  X 

20 

5,635509  +  1 . 5321 IX  <  .15001044X+2  -  .32412528X^3  ■*-  Y 

21 

Y  t-   .04883571X+4  +  .4137684X+5  '<-    .25185536X^6  ■*■   Y 

22 

Y  f  052733125Xt7  +  .002576072X+8  •>  Y;  lOtY  ->  R  (Rib  +  4) 

O"' 

FLT  4;  PRT  RR16;  SPC  1;  PRT  R  ( R 1 6  +  1),  R  (R16  +  2);  SPC  1:  3  ->  X 
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In  this  section  v(t  )  for  both  channels  is  measured.  These  values 
are  then  put  into  the  curve   fit  equations  for  v(t  )  vs  R/hr;  equations  4- 
1  and  2. 

24:  PRT  R  (R16  +  X);  X  +  1  -*■  X,  IF  X  <  4;  JMIJ  0 
!5:  S^C  1;  FXD  4:  PRT  (R  (Rib  +  3)  -R  (R16  +  4))/R  (R16  +  3) 

RR16  ->  X;  FXD  4;  SPC  1;  PRT  (R  (Rib  +  3}  -X)/X,  (R  (R 16  +  4)  -X)/X 

R16  +  5  -»  R16,  SPC  3 

GTO  "A:' 

"1";  R16  ->  RRI7 

ENT  "TAPE  ?",  A,"  FILE  ?\  A 

RCF  A,  R  (R16  -  1) 

GTO  "END" 

"DLL AY" 

The  compensation  values  tabulated  in  tables  4-1  and  -2  are  then  cal- 
culated and  stored. 
34:  CMD  !,ZV5"',  FMT  *;  RED  13,  A 
35:   "0";  CMD  "?V5";  RED  13,  B 
36:   IF  X  f   35:  10  ->  X 
37:   IF  E-A  <  X;  GTO  "0" 
38:  RET 
39:  "EMD":  END 

This  is  a  subroutine  used  to  delay  the  computer  35  see  between  mea- 
surements . 

Plotting  Routine 

The  foil  owing  in  the  data  output  plotting  routine  used  to  produce 
figures  4-8  through  4-23, 


9S 


0:  PRT  "FROG.  PLOTS  I  VS.;' 
1:  PRT  "PIC  VOLTAGE  SIG." 

PRT  "DATA  FROM  FILE" 

PRT  "GIVEN.";  SPG  3 

FNT  "TAPE  ?\  B;  ENT  "PLOT  FILE  ?",  A;  LDF  A 

SGL  -8200,  11 000..  -2000,  11000 

LTR  2000,  10600,  321 

PLT  "CURRENT  VS  PIC  SIG." 
8:  0  ■*■  X  ;  0  ■>  Y 
9:  PEN;  PL!  X,  Y;  1  ->   A;  R9  ■*  B 

1000/ (-LOG  R9  +  LOG  RIO)  ->  Rll 

{•-LOG  R9  +  LOG  (B  *  A))  *  Rll  -  X 

PIT  X,  Y;  PLT  X,  Y  +  100;  PLT  X,  Y 

A  +  I  +  A;  IF  A  <  9;  JMP  -2 

B  *  10  ■>  B;  1  •>  A;  IF  10000  >  (-LOG  R9  +  LOG  (B*A))  *  Rll;  JMP  -3 

PLT  1 0000,  Y 

PEN;  I  -y  A;  R6  ->  B:  0  ->  X;  0  ■+  Y;  PLT  X,  Y 

10000/ (-LOG  R6  +  LOG  R7)  ->  R8 

(-LOG  R6  h  LOG  (B  *  A))  *  R8  ->  Y 

PLT  X,  Y;  PLT  X  +  10  0,  Y;  PLT  Xs  Y 
20:  A  +  1  -*  A;  IF  A  <  9;  JMP  -2 

21:  B  *  10  ->  B;  i  •+  A;  IF  10000  >  (-LOG  R6  +  LOG  (B  *  A))  *  R8;  JMP  -3 
22:  Pi  J  x,  lOQG'O 

10000  *  X;  0  ■*■   Y 

24:  PEN;  PLT  X,  Y;  PLT  X,  Y  *  10000;  PLT  X  -  10000,  Y  +  10000;  PEN 
25:   R9  >  B;  0  -+  Y 
26:   (-LOG  R9  +  LOG  B)  *  Rll  -*  X;  B  *  10  -  3;  IF  B  >  RIO;  JMP  2 


i  ■:. : 

13: 

14; 

15 

16: 

1  7 
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27:  PLT  X,  Y;  PLT  X,  Y  +  10000;  FEN;  JMP  -i 

28:  RC  ->  5;  0  ->  X 

29:  (--LOG  R6  +  LOG  B)  *  R3  ->   Y;  B  *  10  -*  B;  IF  B  >  R7;  JMP  2 

30:  PLT  X,  Y;  PLT  X  +  10000,  Y;  PEN;  JMP  -  ! 

31 :  R6  ->■   B 

32:  -1000  ->   X;  LOG  B  -»■  A 

33:  (--LOG  R6  +  LOG  B)  *  R8  -  100  ->   Y 

34:  LTR  X,  Y,  211:  PLT  "10" 

35:  LTR  X  +  300,  Y  +  300;  FXD  0;  PLT  A;  A  +  1  ->■  A;  IF  A  <  LOG  R7;  10  *  B 

-  3;  JMP  -2 

36:  LTR  3000,  -  900,  211 

37:  PLT  "IONIZATION  CURRENT  (AMPS)" 

38:  R9  -►  B,  -500  ->  Y;  LOG  8  ■?     A 

33:  (-LOG  R9  +  LOG  B)  *  Rll  -300  ■*■   X 

40:  LTR  Xs  Y ,  211;  PLT  "10" 

41:  LTR  X  +  300,  Y  +  300;  FXD  0;  Pi T  A;  A  +  1  ->  A;  IF  A  <  LOG  RIO;  10 

*  B  ->  B;  JMP  -2 

42:  LTR  -1200,  4300,  212 

43.  PLT  "PIC  S1G.  (VOLTS)" 

44:  LTR  -  7300,  10800,  211 

45:  PIT  "IG-SS  VG(+)  IF-SS  VF  (  )" 

46:  -2200  ->  X;  10390  ■+  Y;  PLT  X,  Y;  PEN 

47:  PLT  X  +  50,  Y  +  50;  PLT  X  +  50,  Y  -  50;  PLT  X  -  50,  Y  -  50;  PLT  X  -  50 : 

Y  -  50;  PLT  X  +  50,  Y  +  50;  PEN 

48:  LTR  -7900,  10600,  111 

49:  PLT  "AMPS  VOLTS  AMPS  VOLTS" 


50; 


51 


53: 

54 


57 

58: 

59: 

60: 

61 

62: 

63: 

64: 

65: 

65: 

67: 

58: 

69: 

70: 

7^  , 

/  I  ■ 

11: 
73: 

It  • 


11000  -»  Y;  -600  -  Z;  PIT  -8200,  Y;  FIT  -8200,  Z;  PLT  -6545,  Z;  PLT 

-6545,  Y;  PLT  -4890,  Y 

PLT  -4890,  Z;  PLT  -3235,  Z;  PLT  -3235 .  Y;  PLT  -8200,  Y;  PLT  -1580,  Y; 

PLT  -1580:.  Z;  PLT  -8200 ,  Z  - 

PEN,  PLT  -8200 ,  10480;  PLT  -15805  10480 

10350  ->  Z;  21  *  A;  0  ->■  R12;  A  ■->-   R12 

"0";  ENT  "CHANGE  PEN?1',  C 

Z  -250  -*  7. 

"T'j  LTR  -8100,  Z,  211;  FLT  2;  PLT  RR13 

LTR  -5460,  Z:  FLT  2;  PLT  R  (R13  +  1) 

LTR  -4790,  Z;  FLT  2;  PLT  R  (R13  +  2) 

LTR  -3135,  Z;  FLT  2;  PLT  R  (R13  +  3^ 

Z  -300  ->  Z;  R13  +  4  ■*   R13;  IF  R13  +  1  <  RR^2;  GTO  "1" 

"2":  A  +  1  -  8;  IF  RA  <  1E-14;;  GTO  "3" 

IF  RA  <  R9;  GTO  "3" 

IF  RA  >  RIO;  GTO  "3" 

IF  R8  <  R6;  GTO  "3" 

IF  Rii  >  R7;  GTO  "3" 

(-LOT  R6  ■»  LOG  RB)  *  R8  ■•.  V;  (-LOG  R9  +  LOG  RA)  *  Rll  -*■  X;  PLT  X  -f  50,  V 

PLT  X-50,  Y;  PEN;  PLT  X,  Y  +  50;  PLT  X,  Y  -  50:  PFN 

"3";  A  +  3  -►  C;  A  +  2  ->  B 

TF  R3  <  IE  -  14;  GTO  "4" 

IP  RC  5_  P5:  GTO  "4" 

IF  RC  >  R7;  GTO  "4" 

IF  RB  <_  R9;  GTO  "4" 

IF  RB  >  RIO;  GTO  "4" 

(-LOG  R6  A  LOG  RC)  *  R8  ->  Y ;  (-LOG  R9  h  lOG  RB)  *  Rll  ■*■   X;  PLT  X,  Y; 

PLN,  PLT  X  r   50,  Y  +  50 


77; 
78: 


PL V  X  +  50,  Y  -  50;   PLT  X  ~  50,   Y  -   50;   PLT  X  -  50,   Y  +  50;   PIT 

X  +  50,   Y  +  50,   PEN 

"4";   A  -;   4  -►  A;    IF  A  -»    1    <  RR12;   FTO   "2,: 

R12   <■  1   ->  R12;   IF-  Rl 2  <_  4,   IF  RR12  >  0;   GTO  "9" 

SIP;   END 
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